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•  analysis  of  AFRL  data  for  possible  correlation  with  measurements  made  by  the  Service  d'Aeronomie  (S.A.)  (see  Section  2), 
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•  an  analysis  of  the  ozone  profile  structure  related  to  air  mass  temperature  and  modeling  of  luminance  fluctuation  using  AFRL  and  S.A.  data 
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We  collected  available  measurements  from  S.A.  for  June  98.  The  date  and  location  measurements  of  the  AFRL  data  did  not  satisfactorily 

match  those  of  the  S.A.  data,  so  it  was  difficult  to  make  a  reliable  comparison. 
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1.  SCIENTIFIC  OBJECTIVE 

The  objective  is  the  validation  of  turbulence  and  propagation  codes  at  middle  altitudes/latitudes  and 
the  proof  of  principle  for  flight  of  U.S.  ozone  sondes  and  thermosondes  on  a  major  French  research  balloon 
called  SPIRALE  (SPectrometre  InfraRouge  Atmospherique  a  Laser  Embarque),  successfully  flown  from 
Gap  on  June  29,  1999  and  June  20,  2001.  SPIRALE,  funded  by  CNES  (Centre  National  d’Etudes  Spatiales), 
Onera  and  INSU  (Institut  National  des  Sciences  de  l’Univers),  was  developed  by  Onera  and  LPCE 
(Laboratoire  de  Physique  et  Chimie  de  l’Evironnement)  to  measure  atmospheric  trace  gases  (CO,  03,  CH4, 
N20,  N02,  NO,  C02,  H20,  HF,  etc.)  from  1 3  to  35  km. 

A  test  campaign  (“Gap  98”)  was  conducted  in  June  1998  by  the  AFRL  team  with  AFRL  balloon- 
borne  ozone  sondes  and  thermosondes  above  the  Gap-Tallard  Airfield,  in  southern  France. 

The  town  of  Gap  is  located  approximately  160  km  due  north  of  Toulon,  France.  The  Gap-Tallard 
Airfield  is  situated  in  the  Alps  of  Southeastern  France  on  a  flat  plateau  600  m  above  sea  level.  Small  1200- 
gram  balloons  were  launched  to  obtain  in  situ  vertical  profiles  of  ozone,  wind  speed,  wind  direction, 
humidity,  temperature  and  temperature  fluctuations. 

Dual  ozone  sonde  [l]/thermosonde  [2]  payloads  were  employed  in  some  of  the  measurements.  This 
campaign  is  the  first  in  which  both  ozone  and  temperature  fluctuation  (turbulence)  data  were  simultaneously 
measured  in  situ  from  the  same  small  balloon  platform  at  altitudes  from  ground  level  to  30  km. 

Several  nights  of  sequential  ozone  altitude  profiles  obtained  from  these  in  situ  measurements  showed 
interesting  patterns  above  and  below  the  tropopause  level  (at  12.5  km)  suggestive  of  gravity  waves. 

The  data  are  unique  in  that  they  give  insights  into  not  only  the  structure,  but  also  the  variability  of 
that  structure  over  one-  to  three-hour  periods.  The  integrated  ozone  concentration  agrees  well  with  the  totals 
and  trends  observed  by  the  TOMS  (Total  Ozone  Mapping  Satellite)  satellite  for  a  nearby  reporting  location 
for  the  same  period  (the  altitude-corrected  balloon  data  are  within  4-5%  of  the  TOMS  data.) 

The  AFRL  proposed  to  Onera  to  join  efforts  to  analyze  these  data. 

The  following  tasks,  presented  in  this  report  were  proposed  to  and  carried  out  by  Onera  [3]: 

•  analysis  of  AFRL  data  for  possible  correlation  with  measurements  made  by  the  Service  d’Aeronomie 
(S.A.)  (see  Section  2), 

•  an  example  of  atmospheric  stability  analysis  in  collaboration  with  Dr.  Dalaudier  from  the  Service 
d’Aeronomie.  This  work  is  presented  in  Section  3, 

•  analysis  of  Gap  Cn2  profile  measurements  and  study  of  the  effects  of  atmospheric  turbulence  on  optical 
propagation  in  the  stratosphere.  This  optical  propagation  study  is  reported  in  Section  4, 

•  an  analysis  of  the  impact  of  ozone  profiles  measured  during  the  Gap  campaign  on  the  spectral  radiance  in 
the  8-12  pm  band. 
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2.  DATA  r.  Louvet) 

2.1.  AFRL  data 

The  data  transfer  process  was  completed  on  April  23,  2001. 

The  AFRL  campaign  was  conducted  in  June  1998.  The  list  of  AFRL  flights  including  the  date,  time 
and  type  of  each  balloon  payload  is  given  in  Figure  1. 

The  data  files  received  (see  Figure  2)  included: 

•  wind  data  from  the  first  flight  to  flight  17  on  June  27,  1998, 

•  simultaneous  ozone  and  Cn2  data  measurements  from  flight  9  (on  June  24)  and  flights  10,  1 1  and 
14  (on  June  26), 

•  ozone  data  from  flights  12, 15  and  16. 

A  sample  of  available  measurements  as  a  function  of  the  date  is  given  in  Figure  3. 

The  file  format  is  explained  in  the  dataform  file  given  by  AFRL  (see  Appendix  1). 

This  file  provides  angular  data  and  wind  data  as  a  function  of  altitude  according  to  meteorological 
convention. 

The  ozone  and  temperature  files  provide  measurements  as  a  function  of  altitude.  The  altitude  is 
calculated  from  the  hydrostatic  equation  using  Vaisala  pressure  data.  The  aspect  and  acquisition  frequency 
of  the  data  were  unusual,  so  considerable  additional  information  was  requested  from  AFRL.  The  two  sets  of 
questions  and  the  answers  are  given  in  Appendix  2. 

2.1.1.  Sets  of  data  graphs 

Two  sets  of  graphs  have  been  plotted. 

The  first  set  is  presented  in  Appendix  3.  It  is  intended  to  assist  the  investigation  of  each  flight. 

The  second  set  of  graphs  (Appendix  4)  is  dedicated  to  a  more  complex  analysis  in  relation  with  a 
stability  study.  Here  we  focused  on  runs  9,  10,  11,  and  14  which  provided  the  simultaneous  03/T/Cn  data 
necessary  for  the  stability  analysis. 

In  particular,  we  calculated  the  Richardson  number  and  presented  an  analysis  of  the  principal 
parameters  in  the  three  altitude  ranges:  [14-18]  km  [18-22  ]  km  and  [22-26]  km. 

We  investigated  the  missing  data  link  which  induced  incorrect  shear  values  (flat  values). 

We  detected  some  missing  link  on  data  wind: 

•  run  9  from  15  to  17  km, 

•  run  10  none, 

•  run  1 1  from  3  to  7  km,  from  18.5  to  19.3  and  19.7  km  and  27  km, 

•  run  14  at  19.3  km  and  18.7  km). 
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2.1.2.  Remarks 

The  best  data  with  no  missing  link  data  seem  to  be  from  run  10.  Since  atmospheric  stability  studies 
require  very  good  data  correlation,  run  10  was  chosen  first  to  observe  stability  processes. 

(For  this  work,  see  Section  3) 

2.2.  Available  measurements  provided  by  the  Service  d’Aeronomie 

Several  contacts  were  established  with  researchers  of  the  Service  d’Aeronomie  in  order  to  identify 
common  scientific  interest  and  possible  correlation  measurements.  A  description  of  Service  d’Aeronomie’s 
activities  is  available  on  web  site:  http://www.aero.jussieu.fr/SA. 

Measurements  were  provided  by  five  instruments  available  in  June  1998:  SAOZ,  ozone  lidar,  wind 
lidar,  Rayleigh  temperature  lidar,  and  ozone  sonde.  A  short  description  of  these  measurements  is  given  in 
Appendix  5  and  references  [4-12]. 

The  list  of  available  S.A.  data  as  a  function  of  the  date  and  time  is  shown  in  Figures  3  and  4. 

Only  SAOZ  flights  were  made  from  the  Gap  launch  base.  All  other  measurements  were  made  from 
the  OHP  (Observatoire  de  Haute  Provence).  The  OHP  is  located  about  100  km  south  of  Gap. 

2.2.1.  Comparison  between  AFRL  and  SA  data 

Looking  at  Figures  1,  3  and  4,  we  can  see  that  it  was  unfortunately  difficult  to  provide  a  reliable 
comparison  between  AFRL  data  and  S.A.  data.  Data  were  from  different  dates,  e.g.  SAOZ,  or  from  different 
locations  (AFRL  from  Gap  and  S.A.  from  OHP). 

A  comparison  between  different  ozone  data  profiles  is  given  in  Section  6. 

3.  ANALYSIS  OF  ATMOSPHERIC  STABILITY  (Dr.  Dalaudier  AND  Dr.  Louvet) 


3.1.  Detection  of  predominant  gravity  waves 

AFRL  surveys  (wind,  O3,  H20)  revealed  the  essential  structures  of  the  local  dynamics,  dominated  by 
gravity  waves.  The  wind  zonal  and  meridian  profiles  clearly  show  the  typical  oscillations  associated  with 
these  waves. 

As  an  example,  during  flight  10  in  the  altitude  range  between  18  and  22  km  (see  Figures  5  and  6),  the 
profiles  of  two  wind  components  show  the  signature  of  a  predominant  wave.  A  diagram  of  3D  wind  data 
clearly  shows  the  helicoidal  variation. 

Remark:  outside  the  18-22  km  range,  the  main  part  of  the  fluctuation,  which  does  not  have  a  clear 
oscillating  character,  undoubtedly  corresponds  to  a  random  combination  of  many  gravity  waves  (“not  a 
single  dominant  wave  but  a  mixture  of  many  waves”). 

To  confirm  this  assumption  of  a  predominant  gravity  wave  in  the  1 8-22  km  range,  we  plotted  the 
gradients  of  the  two  available  tracers  (passive  parameters  with  respect  to  dynamics  in  this  altitude  range),  i.e. 
the  potential  temperature  and  the  mixing  ratio  of  ozone. 
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In  fact,  the  vertical  gradient  of  logarithmic  potential  temperature  is  proportional  to  the  Brunt-Vaisala 
frequency. 

Figure  7  shows  that  the  correlation  between  the  ozone  gradient  (grad03)  and  the  Brunt-Vaisala 
pulsation  (N2)  and  the  wind  is  excellent.  Excellent  but  not  perfect  (see  18  km)  because  the  03  gradient  is 
locally  variable  whereas  the  potential  temperature  variations  are  always  positive.  The  same  characteristic 
periods  but  out  of  phase  can  be  observed  on  two  wind  components  and  the  tracers.  The  characteristic 
associated  wavelengths  are  kilometric  in  the  lower  stratosphere. 

3.2.  Study  of  local  stability  conditions 

The  fluctuations  induced  by  the  waves  determine  areas  of  more  or  less  great  stability. 

A  distinction  can  be  made  between: 

Static  stability  (with  respect  to  convective  phenomena)  according  to  a  positive  Brunt-Vaisala 

pulsation  (N2>0),  and  . 

Dynamic  stability,  characterized  by  the  Richardson  number  (Ri=N2/shear  ),  occurring  when  Ri>0.25. 
This  does  not  eliminate  the  possible  presence  of  decreasing  mode  turbulence  in  stable  zones  (see  below),  or 
in  the  opposite  case  of  areas  becoming  dynamically  unstable  where  turbulence  does  not  have  time  to 
develop. 

A  few  interesting  typical  cases  were  observed: 

a)  Possible  contribution  of  atmospheric  layer  (sheets)  [13]. 

Flight  10  at  19.2  km  (see  Appendix  4,  Figure  8-2,  Run  10),  we  observed  a  Cn2  maximum  correlated 
with  a  shear  minimum  and  an  N2  maximum.  This  could  be  the  contribution  of  a  sheet  [1j]  to  the 
temperature  probe  signal,  but  the  origin  of  the  layer  cannot  be  determined  from  the  data. 

b)  Turbulence  associated  with  a  dynamic  instability 

Flight  11  at  an  altitude  of  around  21  km  (see  Appendix  4,  Figure  8-2,  Run  11),  we  observed  a  Cn 
peak,  with  a  high  value  of  N2  and  a  shear  maximum  which  could  be  a  dynamic  instability 

c)  Fluctuation  in  a  neutralized  layer 

Flight  9  at  around  18.1  km  (see  Appendix  4,  Figure  8-2,  Run  9),  a  negative  or  weak  N2  was 
associated  with  fluctuation  of  cn2.  Two  explanations  are  possible: 

The  existence  of  fluctuation  in  a  stable  layer  neutralized  by  a  mixing  process. 

A  possible  artefact  (see  [14])  or  fluctuation  on  a  small  scale  not  observed  by  the  main  Vaisala 
temperature  sensor. 
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3.3.  Conclusion  of  the  atmospheric  stability  study 

Sample  cases  of  analysis  and  interpretation  with  emphasis  on  the  dynamical  state  of  the  atmosphere 
were  presented.  The  general  domination  of  large  scale  dynamics  by  gravity  wave  was  confirmed,  and  a  short 
case  study  of  a  predominant  monochromatic  gravity  waves  was  presented.  Various  typical  atmospheric 
stability  cases  were  observed  within  this  short  example  of  stability  analysis.  The  numerous  peaks  observed 
with  the  thermosonde  are  not  necessarily  associated  with  atmospheric  turbulence  [14], 

3.4.  Prospects 

This  work  is  a  preliminary  example  of  possible  investigation.  The  following  future  work  is 
proposed. 

A  more  detailed  and  more  complete  analysis  of  the  existing  layers  (sheets),  with  investigation  of  their 
origin  and  their  influence  on  the  Cn2,  could  be  carried  out. 

An  analysis  of  the  horizontal  structure  and  its  variation  over  time  could  be  made  by  correlation 
between  the  closed  flights  as  well  as  an  investigation  of  the  variation  of  gravity  waves  and  turbulent  layers. 

A  statistical  stability  analysis  of  the  static  and  dynamic  layer  with  a  high  Cn2  level  would  make  it 
possible  to  better  identify  the  real  contribution  of  turbulence  on  the  signal. 


4.  OPTICAL  PROPAGATION  (Dr.  Michau  and  Rousset-Rouviere) 

This  chapter  is  aimed  at  describing  the  effects  of  atmospheric  turbulence  on  the  propagation  of 
optical  beams  over  a  long  path  (500  km)  in  the  upper  atmosphere  (1 0  to  25  km). 

This  study  was  made  by  conducting  numerical  experiments  under  various  conditions  of  propagation. 
The  first  part  of  this  section  describes  the  model  used  for  the  numerical  experiments.  The  second  part 
presents  the  conditions  of  the  different  experiments.  The  results  of  the  numerical  experiments  are  presented 
in  the  last  part  of  this  section. 

4.1.  Description  of  wave  propagation  through  atmospheric  turbulence 

4.1.1.  Plane  wave  propagation 

The  model  describing  propagation  of  a  plane  wave  through  atmospheric  turbulence  is  based  on  a 
method  widely  described  in  the  literature  [15].  In  this  method,  the  phase  delays  induced  by  the  spatial 
variations  of  the  refraction  index  are  described  using  phase  screens  placed  regularly  along  the  propagation 
axis.  The  statistical  properties  of  the  optical  delay  map  of  each  phase  screen  has  a  stationary  Gaussian 
distribution,  deduced  from  Kolmogorov’s  law.  The  electromagnetic  field  EMF  is  propagated  along  the 
distance  between  each  phase  screen  using  the  propagation  equation  in  a  vacuum,  i.e.  with  a  Fresnel 
transform. 

For  a  plane  wave,  the  propagation  between  phase  screen  q  and  phase  screen  q+1  is  written  (see 
Figure  8): 
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Uq+i  (©)=  t/'+1  (a))exp(-i7tAh6)2 ).exp (ikh) 

where  Uq (&>)  is  the  Fourier  transform  of  the  EMF  complex  amplitude  U q[p)  in  plane  q  before  phase 
screen  q  and  U'q{p)  is  the  EMF  complex  amplitude  after  phase  screen  q,  denoted  (pturhq : 

U'q{p)=Uq{p)exp  ^(pturbXp} 


h  denotes  the  distance  between  phase  screens  q  and  q+1,  X  the  wavelength  and  k  -  — . 

4.1.2.  Spherical  wave  propagation 

The  model  presented  in  Section  4.1.1  was  modified  to  simulate  spherical  wave  propagation.  The 
EMF  complex  amplitude  is  then  broken  down  as: 

u(p)=m(p)usw{p) 

where  f/w(p)  represents  the  EMF  complex  amplitude  corresponding  to  a  perfect  spherical  wave  and  m  a 
modulation  term  taking  into  account  the  perturbations  due  to  turbulence  effects.  With  this  definition,  the 
modulation  term  obeys  the  following  propagation  equation: 

mq+ 1  (fl>)=  fnq  (®)exp(-  ikXco2heq ) 

where  h  =  — - — ,  h  is  the  distance  between  phase  screens  q  and  q+1,  and  z  is  the  distance  from  the 
eq  \-h/z 

source  to  phase  screen  q+1.  This  allows  propagation  of  a  spherical  wave  to  be  described  in  the  same  way  as 
propagation  of  a  plane  wave.  The  beam  geometry  is  shown  in  Figure  9.  Note  that  the  sizes  of  the  phase 
screens  increase  linearly  with  the  propagation  distance.  Turbulence  effects  between  the  source  and  the  first 
phase  screen  are  not  taken  into  account  in  our  model. 

Using  these  models,  a  software  code  was  developed  at  Onera  to  describe  propagation  of  plane  and 
spherical  waves  through  atmospheric  turbulence.  The  phase  screens  are  drawn  randomly.  The  main  input 
parameters  of  this  code  are: 

-  the  conditions  defining  turbulence  along  the  optical  path:  the  Cn2  profile,  the  outer  scale  (Lo),  the  inner 
scale  (lo), 

-  the  beam  characteristics:  wavelength,  distance  between  the  source  and  the  receiver  plane,  width  of  the 
receiver  plane, 

-  the  numerical  parameters:  spatial  sampling  of  the  phase  screens,  number  of  pixels  in  a  screen  phase, 
position  of  the  phase  screens  along  the  propagation  axis. 
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The  first  output  of  the  software  is  the  EMF  complex  amplitude  in  the  receiver  plane.  For  a  wave 
propagated  through  turbulence,  the  EMF  complex  amplitude  in  the  receiver  plane  can  be  written: 

E(r)  =  E0(r)exp[iy/(r)] 

E0  is  the  EMF  complex  amplitude  without  perturbations  and  exp[iy/(p)]  represents  the  perturbation  term 
due  to  turbulence  effects,  y/  is  a  complex  number  whose  real  part  represents  the  phase  perturbations  denoted 
(p  and  whose  imaginary  part,  the  so-called  log-amplitude,  represents  the  amplitude  perturbations,  i.e.  the 
scintillation  effects: 


Kr)  =  <p(r)  +  iz(r) 

A  numerical  experiment  was  conducted  by  processing  a  set  of  EMF  complex  amplitudes,  Ej,  obtained 
in  the  receiver  plane  with  the  same  set  of  input  parameters  but  with  different  drawings  of  the  phase  screens. 
Then,  the  statistical  properties  of  (p  and  y/axz  determined  from  the  set  of  Ej  values: 

-  the  spatial  power  spectral  density  (PSD)  of  cp , 

-  the  spatial  variance  of  y/  as  a  function  of  the  propagation  distance, 

-  the  spatial  PSD  of  y/ , 

-  the  spatial  PSD  of  I=|E|2. 

4.2.  Propagation  conditions 

This  part  describes  the  different  propagation  conditions  that  were  simulated  for  the  numerical 
experiments.  The  propagation  geometry  is  shown  in  Figure  10: 

-  the  receiver  plane  was  1 5  km  high, 

-  the  source  height  varied  between  10  and  25  km, 

-  the  propagation  distance  was  500  km. 

The  outer  scale  was  assumed  to  be  50  meters  and  the  inner  scale  2  mm.  We  studied  propagation  at 
the  following  wavelengths:  1.06,  1.55,  5  and  10  pm.  Turbulence  was  represented  by  20  phase  screens 
sampled  with  a  256  by  256  grid.  The  size  of  the  last  screen  was  50  meters. 

The  Cn2  profile  used  in  the  numerical  experiment  was  obtained  by  a  log-linear  fit  of  the  data  given  in 
the  GAP98 1 1  .txt  file.  The  equation  of  fitted  profile  is: 

C£  =cr.l0H*) 

where  h  is  the  altitude  and  a  and  (i  are  constants.  The  fit  was  made  in  the  [5;  25  ]  km  range. 

The  figure  1 1  shows  the  C„2  profile  corresponding  to  the  GAP98 1 1  .txt  file  and  the  result  of  the  fit. 
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4.3.  Results 

Table  1  shows  the  main  parameters  characterizing  optical  propagation  for  the  numerical  experiments: 

-  the  geometry  of  the  optical  path, 

-  the  Fresnel  distance  ( *J~XL ), 

-  the  Fried  parameter  ro, 

-  the  standard  deviation  of  log-amplitude  %  in  the  Rytov  regime  (weak  perturbations), 

-  the  standard  deviation  of  log-amplitude  %  estimated  with  the  numerical  experiment. 

As  indicated  in  table  1,  the  log-amplitude  standard  deviation  <JX  obtained  from  the  Rytov 
approximation  and  the  one  obtained  from  the  numerical  experiment  are  very  similar,  except  for  the  worst 
case  i.e.  for  hsource  =  10  km,  X  =  1.06  |im.  In  this  case,  the  standard  deviation  ax  obtained  from  the  numerical 
experiment  is  smaller  than  the  Rytov  one,  c  .  This  result  indicates  the  beginning  of  the  saturation 

regime  of  scintillation.  In  the  other  cases,  the  numerical  experiments  were  conducted  in  the  weak 
perturbation  regime.  This  behavior  might  be  able  to  be  predicted  by  comparing  r0  with  the  Fresnel  distance 

4al. 


Table  1 


Source  height  (km)  -  Wavelength 
(l-im) 

s[XL  (m) 

ro(m) 

ax 

^  X  Rylov 

10-1.06 

0.73 

0.39 

0.74 

0.84 

10-1.55 

0.88 

0.61 

0.68 

0.67 

10-5.00 

1.58 

2.51 

0.35 

0.34 

10-10.0 

2.23 

5.76 

0.22 

0.22 

15-1.06 

0.73 

0.47 

0.65 

0.63 

15-1.55 

0.88 

0.74 

0.54 

0.51 

15-5.00 

1.58 

3.05 

0.26 

0.25 

15-10.0 

2.23 

7.0 

0.17 

0.17 

20-1.06 

0.73 

0.56 

0.47 

0.45 

20-1.55 

0.88 

0.88 

0.41 

0.40 

20  -  5.00 

1.58 

3.60 

0.20 

0.20 

20-10.0 

2.23 

8.26 

0.13 

0.13 

25-1.06 

0.73 

0.64 

0.39 

0.39 

F  25-1.55 

0.88 

1.01 

0.31 

0.31 

25  -  5.00 

1.58 

4.13 

0.15 

0.16 

25  -  10.0 

2.23 

9.5 

0.10 

0.10 

The  results  obtained  from  the  numerical  experiments  are  given  in  Figure  12  to  Figure  27.  Each  figure 

shows: 

-  an  example  of  intensity  distribution  in  the  receiver  plane, 
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-  the  spatial  power  spectral  density  (PSD)  of  (p , 

-  the  spatial  standard  deviation  %  as  a  function  of  the  propagation  distance  or  as  a  function  of  the  value 
obtained  from  the  Rytov  approximation, 

-  the  spatial  PSD  of  %, 

-  the  spatial  PSD  of  I=IE12. 

For  the  weak  perturbation  case,  the  PSD  of  %  and  I  and  the  standard  deviations  obtained  from  the 
numerical  experiments  were  fitted  with  analytical  results  obtained  from  the  Rytov  model.  Saturation  of  the 
intensity  fluctuations  can  be  seen  in  Figure  12.  The  PSD  of  I  obtained  from  the  numerical  experiment  is 
lower  than  the  Rytov  one  mainly  for  intermediate  spatial  frequencies  [16,  17]. 

The  examples  of  intensity  distribution  in  the  pupil  plane  of  the  receiver  also  illustrate  this  behavior. 
A  single  characteristic  structure  size  (the  Fresnel  distance)  can  be  seen  in  the  weak  perturbation  case.  In  the 
saturation  regime  (Figure  12)  filaments  appear  and  two  structure  sizes  can  be  distinguished:  the  width  and 
the  length  of  these  filaments. 

4.4.  Conclusion 

Numerical  experiments  were  conducted  to  simulate  the  effects  of  atmospheric  turbulence  on  the 
propagation  of  optical  beams  over  a  long  path  (500  km)  in  the  upper  atmosphere  (10  to  25  km).  The  study 
showed  that  propagation  could  be  analyzed  using  the  weak  perturbation  model  (Rytov  approximation)  in 
most  cases.  In  the  worst  case,  (a,=  1.06  |im,  hsourCe=  10  km,  hreceiver  =  15  km,  distanceSOurce-receiver  =  500  km), 
the  intensity  fluctuations  entered  the  saturation  regime.  The  PSD  of  the  phase  and  log-amplitude  were 
characterized. 

This  study  was  made  using  a  phase  perturbation  model  based  on  the  perturbation  of  the  refraction 
index  in  the  inertial  range.  This  model  is  probably  valid  for  describing  the  intensity  variations  inside  the 
receiver  pupil.  However,  it  does  not  take  into  account  large  scale  fluctuations  of  the  refraction  index  along 
the  optical  path  such  as  gravity  waves.  These  variations  can  substantially  modify  the  amplitude  of  the 
variations  of  the  total  intensity  in  the  receiver  pupil. 
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5.  ANALYSIS  OF  THE  OZONE  PROFILE  AND  RADIANCE  COMPUTATIONS 

(Dr.  Simoneau) 

The  goal  of  this  study  is  to  analyze  the  impact  of  ozone  profiles  measured  during  the  Gap  campaign 
on  the  spectral  radiance  in  the  8-12  pm  band.  In  this  work,  we  run  the  radiative  transfer  code  MATISSE 
developed  at  Onera,  using  ozone,  temperature  and  water  vapor  profiles  measured  during  the  campaign  as 
input  data. 

The  5.1.  paragraph  contains  a  brief  description  of  MATISSE.  It  is  followed  by  a  discussion  on 
profiles  measurements  in  section  5.2.  Computation  results  are  finally  discussed  in  paragraph  5.3. 

5.1.  MATISSE  code 

MATISSE  1.1  (MATISSE  stands  for  “Advanced  Earth  Modeling  for  Imaging  and  Scene 
Simulation”)  is  the  first  version  of  a  new  program  currently  under  development  at  Onera  [18],  whose 
purpose  is  to  compute  background  radiance  images  by  taking  into  account  atmospheric,  clouds  and  ground 
radiation  as  well  as  the  variability  of  atmospheric  properties  along  the  lines  of  sight.  The  development  is 
planned  for  the  code  to  be  delivered  by  the  mid  year  2002.  Before  MATISSE  1.1  development  phase,  a 
prototype  has  been  realized  whose  first  objective  was  to  compute  transmission  and  atmospheric  radiance 
along  an  optical  path  in  an  atmosphere  exhibiting  geographically  variable  thermodynamic  properties  with  a 
fast  computing  method  for  the  radiation  propagation  (band  model  with  1  cm'1  spectral  resolution).  This 
prototype  has  been  used  for  this  study.  Accordingly,  it  allows  computations  in  atmospheres  defined  by  a 
large  number  of  atmospheric  layers,  which  is  not  the  case  for  most  of  the  available  radiative  transfer  codes. 
We  use  MATISSE  in  the  8-12  pm  spectral  band  without  taking  into  account  aerosols,  and  by  computing 
atmospheric  thermal  emission  only.  Refraction  effects  are  also  taken  into  account. 

5.2.  The  atmospheric  profiles 

Figure  28  and  Figure  29  show  the  ozone  profiles  versus  altitude  measured  during  gajo  campaign.  Four 
ozone  profiles  measured  by  the  Air  Force  Research  Laboratory  (AFRL)  during  the  24th,  25*  and  26th  of  June 
1998  and  the  MidLatitude  Summer  profile  (MLS)  [19]  are  plotted  on  the  Figure  28.  As  can  be  seen,  the 
measured  profiles  are  close  to  the  MLS  profile,  though  exhibiting  small  variations  in  the  vicinity  of  15  km. 
Figure  29  shows  ozone  profiles  measured  by  the  ozone  Lidar  (22nd,  23rd,  25th,  26th  and  29th  of  June),  the  ECC 
probe  (24th  of  June)  and  SAOZ  (23th  of  June)  instruments  (see  Section  2  and  Figure  3,  measurements  of 
Service  d’Aeronomie).  Huge  variations  of  the  Lidar  measurements  are  observed  between  40  and  50  km,  so 
the  major  part  of  this  work  is  devoted  to  the  radiance  variations  induced  by  the  use  of  these  atmospheric 
ozone  profiles  as  input  data.  The  striking  evolution  of  the  SAOZ  profile  below  15  km  is  a  modeling  artifact 
due  to  the  absence  of  measurement  for  these  altitudes. 

Moreover,  the  profiles  are  only  available  for  a  restricted  section  of  the  atmosphere,  so  they  have  been 
extrapolated  from  the  ground  to  the  top  of  the  atmosphere,  and  the  missing  atmospheric  molecules  have  been 
added  by  using  the  PRFL  code  [20].  The  deduced  profiles  are  plotted  on  the  Figure  30  as  well  as  the  MLS 
profile  and  the  climatological  profile  used  in  the  PRFL  code  (so  called  Climato  20/06  on  the  figure).  It  can 
be  seen  that  the  large  fluctuations  of  the  Lidar  profiles  induce  artifacts  in  the  extrapolated  data  around 
50  km. 
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5.3.  Results 

The  figure  31  shows  radiance  computations  performed  by  using  ozone  Lidar  measurements  and  the 
MLS  profile  as  input  data.  The  following  observational  conditions  are  used :  the  observer,  located  at  an 
altitude  of  50  km  is  looking  towards  the  ground  and  the  path  length  is  10  km.  The  figure  put  in  evidence  a 
high  variability  of  the  radiance.  The  largest  gap  between  Lidar  measurements  and  MLS  profile  corresponds 
to  the  data  of  the  23rd  of  June,  where  the  value  at  50  km  is  close  to  zero.  Nevertheless,  this  altitude  is  the 
upper  limit  of  the  Lidar  range,  so  this  result  has  to  be  taken  with  care. 

Figure  32  shows  radiance  computations  using  the  reverse  geometry  :  the  observer  is  now  located  at 
an  altitude  of  40  km  and  is  looking  upward.  The  same  conclusions  as  previously  can  be  drawn,  but  it  appears 
a  less  pronounced  radiance  variability. 

Figure  33  to  Figure  36  show  radiance  computations  for  an  horizontal  line  of  sight,  a  path  length  of 
1  km  and  an  observer  located  at  different  altitudes  on  each  figure.  The  data  used  for  the  computations  are 
again  the  Lidar  profiles  measurements.  The  aim  of  these  computations  is  to  study  the  radiance  variations 
induced  by  using  various  atmospheric  profiles  as  they  could  be  measured  by  an  observer  in  the  atmosphere 
looking  at  a  target  flying  at  the  same  altitude.  For  Figure  33  the  observer  is  located  at  an  altitude  of  50  km. 
As  expected,  the  lowest  radiance  value  is  for  measurements  performed  on  the  23rd  of  June  where  the  ozone 
concentration  measurement  is  close  to  zero.  It  can  be  seen  that  in  all  cases,  the  MLS  profile  induce  the 
highest  radiance  value.  This  result  could  be  unexpected  regarding  the  ozone  value  for  this  altitude  in 
comparison  with  the  other  Lidar  measurements,  but  it  appears  that  the  temperature  value  in  the  MLS  case  is 
higher  than  the  data  recorded  during  the  campaign.  This  last  point  explains  the  higher  radiance  value  for 
these  wavelengths. 

For  the  case  treated  in  the  Figure  34,  the  observer  is  located  at  an  altitude  of  48  km,  which  is  the 
altitude  of  the  ozone  peak  value  for  the  22nd  of  June.  As  expected,  the  radiance  reaches  its  higher  value  for 
this  day.  For  the  other  days,  the  MLS  profile  always  induces  the  higher  radiance  value  for  the  same  reasons 
as  previously  discussed. 

Figure  35  exhibits  a  less  pronounced  variability  in  correlation  with  the  ozone  profile  fluctuations  for 
this  altitude,  and  has  a  low  value  at  an  altitude  of  40  km  where  the  Lidar  measurements  are  smoother 
(Figure  36). 

Figure  37  shows  radiance  computation  using  similar  observation  geometry  as  for  the  Figure  35,  but 
the  path  length  is  100  km.  In  this  case,  the  ozone  concentrations  are  averaged  over  the  optical  path  due  to  the 
sphericity  of  the  atmospheric  layers  inducing,  as  can  be  observed,  a  lower  variability  for  the  computed 
radiance. 

The  results  plotted  on  the  Figure  38  are  obtained  with  another  geometrical  conditions:  in  this  case, 
the  observer  is  located  at  an  altitude  of  46  km  (in  the  validity  domain  of  the  Lidar  measurements)  looking  at 
the  ground  with  a  zenithal  angle  of  135°.  The  final  altitude  of  the  line  of  sight  is  3.7  km.  The  SAOZ  artifact 
already  mentioned,  induces  the  high  radiance  value  in  the  ozone  emission  band.  Nevertheless,  a  line 
structure  of  radiance  values  appears,  which  is  due  to  the  presence  of  water  vapor  at  low  altitude.  This 
phenomenon  is  well  marked  for  the  AFRL  measurement  as  it  can  be  expected,  because  this  profile 
(measured  at  the  same  time)  is  the  most  humid  as  can  be  seen  on  the  Figure  39. 

The  last  two  figures  (Figure  40  and  Figure  41)  show  a  comparison  of  the  computed  radiance  from 
measurements  performed  in  the  same  period  and  with  two  different  devices  (Lidar/SAOZ  and  AFRL/ECC 
sounds).  It  appears  a  marked  gap  between  Lidar  and  SAOZ  measurements,  nevertheless  it  is  not  possible  to 
conclude  because  measurements  have  been  performed  on  two  sites  separated  by  100  km. 
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5.4.  Conclusion 

In  this  study  the  ozone  profiles  measured  during  the  Gap  campaign  which  hold  in  June  1998  have 
been  used  for  estimating  their  impact  on  atmospheric  radiance  computations.  Some  important  variations  can 
be  observed  for  particular  geometries,  when  using  the  high  fluctuating  Lidar  profiles  measurements. 
Nevertheless,  it  is  difficult  to  compare  radiance  variations  induced  by  the  choice  of  the  profile  measurement 
device  because  the  available  data  have  not  been  recorded  on  the  same  site  location  or  at  the  same  time. 


6.  CONCLUSION 

We  collected  available  measurements  from  S.A.  for  June  98.  The  date  and  location  measurements  of 
the  AFRL  data  did  not  satisfactorily  match  those  of  the  S.A.  data,  so  it  was  difficult  to  make  a  reliable 
comparison. 

The  analysis  of  AFRL  data  established  the  existence  of  predominant  monochromatic  gravity  waves 
and  possible  cases  of  dynamic  effects  and  sheets  have  been  observed.  Suggestions  for  future  work  are  given. 

Numerical  experiments  were  conducted  to  simulate  the  effects  of  atmospheric  turbulence  on  the 
propagation  of  optical  beams  over  a  long  path  (500  km)  in  the  upper  atmosphere  (10  to  25  km).  The  study 
showed  that  the  propagation  could  be  analyzed  using  the  weak  perturbation  model  (Rytov  approximation)  in 
most  cases  Future  work  can  be  proposed  to  take  into  account  the  influence  of  refraction  index  fluctuations 
caused,  for  instance,  by  gravity  waves. 

The  Gap98  ozone  profiles  measured  have  been  used  for  estimating  their  impact  on  atmospheric 
radiance  computations.  Nevertheless,  it  is  difficult  to  compare  radiance  variations  induced  by  the  choice  of 
the  profile  measurement  device  because  the  available  data  have  not  been  recorded  on  the  same  site  location 
or  at  the  same  time. 
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Figure  1  -  List  of  flights  from  the  GAP98  campaign  provided  by  AFRL. 
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Figure  2  -  Available  AFRL  data  files. 
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Figure  3  -  Available  data  from  AFRL  and  S.A.  as  a  function  of  date. 
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Figure  8  -  Diagram  of  plane  wave  propagation  along  z  between  phase  screens  q  and  q+1.  U’q  is  the  EMF  complex  amplitude 
after  phase  screen  q  and  Ul/+1  is  the  complex  amplitude  before  phase  screen  q+1. 
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Figure  9  -  Diagram  of  spherical  wave  propagation  along  z ■  Zq  is  distance  between  the  first  phase  screen  and  the  source 
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Figure  10  -  Propagation  geometry. 
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Figure  12  -  Result  of  the  numerical  experiment  (hsource  =  10  km,  A. =1.06  pm) 
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Figure  13  -  Result  of  the  numerical  experiment  (hwurce  =  10  km,  X  =  1.55  pm). 
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Figure  14  -  Result  of  the  numerical  experiment  (hsource  =  10  km,  A, -5.0  pm). 
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Figure  15  -  Result  of  the  numerical  experiment  (hS(mrce  =  10  km,  X=10.0  pm). 

Upper  left  :  PSD  of  <p  (phase  of  the  EMF  complex  amplitude) 

Upper  middle :  PSD  of  X  (log-amplitude  of  the  EMF  complex  amplitude) 

Upper  right  :  PSD  of  I  (square  modulus  of  the  EMF  complex  amplitude) 

Middle  left  :  Variance  of  X  as  function  of  the  propagation  distance 

Middle  right  :  Variance  ofx  as  function  of  the  theoretical  variance  of  X  (Rytov 

approximation ) 

Bottom  :  Example  of  intensity  distribution  in  the  receiver  plane 
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Figure  17  -  Result  of  the  numerical  experiment  ( h  wurce  -  15  km,  X  =  1.55  pm). 


Upper  left 
Upper  middle  : 
Upper  right  : 
Middle  left 
Middle  righ  t: 

Bottom 


PSD  of  cp  (phase  of  the  EMF  complex  amplitude) 

PSD  of  x  (log-amplitude  of  the  EMF  complex  amplitude) 

PSD  of  I  (square  modulus  of  the  EMF  complex  amplitude) 
Variance  ofx  as  function  of  the  propagation  distance 
Variance  ofx  as  function  of  the  theoretical  variance  of  X  (Rytov 
approximation) 

Example  of  intensity  distribution  in  the  receiver  plane 


ONERA 


RTS  2/05199  DOTA 


-38 


SANS  MENTION 
DE  PROTECTION 

SEPTEMBRE  2001 

SMP 


Figure  18  -  Result  of  the  numerical  experiment  (hmurce  =  15  km,  X  =5.0  pm). 
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Figure  19:  Result  of  the  numerical  experiment  (hsource  =  15  km,  A,  =  10.0  pm). 
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Figure  20  -  Result  of  the  numerical  experiment  (hmurce  =  20  km,  X  =  1.06  pm). 
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Figure  21  -  Result  of  the  numerical  experiment  (hwurce  =  15  km,  X  =  1.55  pm). 
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Figure  22:  Result  of  the  numerical  experiment  (hwurce  =  15  km,  X  =5.0  pm). 

Upper  left  :  PSD  of  rp  ( phase  of  the  EMF  complex  amplitude) 

Upper  middle :  PSD  of  X  (log-amplitude  of  the  EMF  complex  amplitude) 

Upper  right  :  PSD  of  I  (square  modulus  of  the  EMF  complex  amplitude) 

Middle  left  :  Variance  of  x  as  function  of  the  propagation  distance 
Middle  right  :  Variance  ofx  as  function  of  the  theoretical  variance  of  X  (Rytov 
approximation) 

Bottom  :  Example  of  intensity  distribution  in  the  receiver  plane 
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Figure  23  -  Result  of  the  numerical  experiment  (hwurce  =  15  km,  X- 10.0  pm). 
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Figure  24  -  Result  of  the  numerical  experiment  (hwurce  =  25  km,  X- 1.06  pm). 
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Figure  25  -  Result  of  the  numerical  experiment  (hsource  =  25  km,  A,  =1.55  pm). 

Upper  lef  t:  PSD  of  (p  ( phase  of  the  EMF  complex  amplitude) 

Upper  middle :  PSD  of  %  (log-amplitude  of  the  EMF  complex  amplitude) 

Upper  right  :  PSD  of  I  (square  modulus  of  the  EMF  complex  amplitude) 

Middle  left  :  Variance  of  X  as  function  of  the  propagation  distance 
Middle  right  :  Variance  of  x  as  function  of  the  theoretical  variance  of  x  (Rytov 

approximation) 

Bottom  :  Example  of  intensity  distribution  in  the  receiver  plane 
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Figure  26  -  Result  of  the  numerical  experiment  (hsoura,  =  25  km,  X  -  5.0  pm). 

Upper  left  :  PSD  of  (p  (phase  of  the  EMF  complex  amplitude) 

Upper  middle :  PSD  of  x  (log-amplitude  of  the  EMF  complex  amplitude ) 

Upper  right  :  PSD  of  I  (square  modulus  of  the  EMF  complex  amplitude) 

Middle  left  :  Variance  of  x  as  function  of  the  propagation  distance 
Middle  right  :  Variance  of  x  as  function  of  the  theoretical  variance  of  X  (Rytov 
approximation) 

Bottom  :  Example  of  intensity  distribution  in  the  receiver  plane 
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obs  =  40km  /  oza=0  /range  =  10km 


wavelength  (cm-1) 


obs  =  50km  /  oza=90  /range  =  1  km 


wavelength  (cm-1) 


obs  =  30km  /  oza=135  /final  altitude  =  12  km 


wavelength  (cm-1) 


obs  =  30km  /  oza=1 35  /final  altitude  =  1 2  km 


wavelength  (cm-1) 
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APPENDIX  1 

FORMAT  OF  AFRL  DATA  FILE 


Gap,  France  1998  Data  to  be  sent  to  France  under  the  Bilateral  US/France  Data  Exchange  Agreement 

Prepared  by  Edmund  A.  Murphy/  AFRL/VSBL,  29  Randolph  Rd.,  Hanscom  AFB,  MA  01731, 
Tel:  (781)377  4403,  e-mail  Edmund. Murphyghanscom.af .mil 
27  March  2001 


NOTE!  -  The  following  information  explains  the  output  formats  for  each  of  two 
balloon  payload  configurations  launched  during  the  Gap,  France  1998  ozone/  thermosonde 
program.  The  flights  consisted  of  either  a  dual  payload  for  the  measurement  of  ozone 
and  turbulence  or  just  the  ozone  experiment.  The  thermosonde  data  is  the  temperature 
difference  across  two  fine  wire  probes  spaced  one  meter  apart.  The  output  data  name  cn2 
is  the  value  of  the  atmospheric  structure  parameter  calculated  from  the  temperature 
difference  across  the  probes.  Ozone  is  given  in  the  output  file  as  Partial  pressure  as 
well  as  number  density. 


There  are  seven  Ozone  payloads : 


1.)  gap9809 

obtained 

on 

06/24/98 

at 

2345 

GMT 

to 

a 

height 

of 

26.6 

km 

is 

a 

dual 

payload . 
2 . )  gap9810 

obtained 

on 

06/25/98 

at 

2130 

GMT 

to 

a 

height 

of 

28.6 

km 

is 

a 

dual 

payload. 
3.)  gap9811 

obtained 

on 

06/25/98 

at 

2359 

GMT 

to 

a 

height 

of 

30.6 

km 

is 

a 

dual 

payload. 

4. )  gap9812  obtained  on  06/26/98  at  1935  GMT  to  a  height  of  30.6  km  is  an  ozone 

payload  only. 

5. )  gap9814  obtained  on  06/26/98  at  2241  GMT  to  a  height  of  30.6  km  is  a  dual 

payload. 

6. )  gap9 815  obtained  on  06/27  98  at  1730  GMT  to  a  height  of  21.6  km  is  an  ozone 

payload  only. 

7. )  gap9  816  obtained  on  06/27/98  at  1926  GMT  to  a  height  of  30.6  km  is  an  ozone 

payload  only. 

8. )  ’ 

The  following  lines  of  code  are  are  taken  from  the  data  reduction  programs  to 
better  explain  the  data  output  files  [gap98xx.txt] 

pro  va_du09g 

;  Ver  oct  98  created  by  EAM  to  process  Vaisala  Thermosondes/ Ozones ondes 

.********************************************************************** 

;  NOTE!!  THIS  PROGRAM  IS  USED  TO  REDUCE  A  DUAL  (THERMOSONDE/OZONE)  PAYLOAD 

;  INPUT  TO  THIS  PROGRAM  IS  OUTPUT  FROM  FULLxxg.PRO 

.********************************************************************** 

/ 

; Create  a  data  .TXT  file 

.  *********************** 

/ 

close, 4 
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openw, 4 , txtf ileout 

for  i=0,fin  do  begin 

printf , 4, format=" (4 (f 12. 4, lx) , 4(el5.4,lx) ,3(el3.5,lx) ) " , zz (i) ,pr (i) , te (i) ,h 
u  (i)  , ecn2 (i) , $ 

parp3 ( i ) , parp4 ( i ) , parp5 ( i ) , numd3 ( i ) , numd4 ( i ) , numd5 ( i ) 

endf or 
close ,  4 

pro  va_ozl2g 

;  Ver  oct  98  created  by  EAM  to  process  Vaisala  Ozonesonde  only 

.************************************************************* 

/ 

;  NOTE! !  THIS  PROGRAM  IS  USED  TO  REDUCE  An  OZONE  PAYLOAD 

;  INPUT  TO  THIS  PROGRAM  IS  OUTPUT  FROM  ozxxg.PRO 

.*************************************** 

t 

;  NOTE!  -  All  five  channels  are  Ozone  * 

,-Calculate  the  Ozone  concentration 
.********************************* 

numdl  =  (7 . 24313*parpl/kte) *1 .e+12 

;This  is  from  Komhyr's  paper  and  is  in  (l/cm/"3)  -  mult,  by  10e+06  for  (l/mA3) 

.  ******************** 

/ 

numdl  =  numdl * 1 . 0  e+ 0  6 

.  ******************** 
t 

;Create  a  data  .TXT  file 
.  *********************** 

close, 4 

openw, 4, txtf ileout 

for  i=0,fin  do  begin 

printf, 4, format=" (4 ( f 12 . 4 , lx) , 6 (el3 . 5 , lx) ) " , zz (i) ,pr ( i ) , te ( i ) , hu ( i ) , $ 
parpl ( i ) , parp2 ( i ) , parp3 ( i ) , parp4 ( i ) , parp5 ( i ) , numdl ( i ) 

endfor 
close, 4 


VARIABLE  NAME 

PARAMETER 

UNITS 

ZZ 

Altitude 

km 

pr 

Pressure 

mbars 

te 

Ambient  temperature 

Deg.  C 

hu 

Relative  humidity 

% 

cn2 

Atmospheric  Structure 

Function 

mA- (2/3) 

parpl 

Ozone  partial  pressure  (channel  #1) 

Nanobars 

numdl 

Ozone  number  density 

(channel  #1) 

/mA3 
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10.  WIND  DATA 


Column  #1  -  Altitude  in  km 

Column  #2  -  Velocity  in  m/sec 

Column  #3  -  Direction  in  degrees  from  true  north 


File  name 

Launch  Date 

Gap9801.win 

06/20/98 

Gap9802.win 

06/21/98 

Gap9804.win 

06/21/98 

Gap9805.win 

06/22/98 

Gap9806.win 

06/22/98 

Gap9807.win 

06/23/98 

Gap9808.win 

06/24/98 

Gap9809.win 

06/24/98 

Gap9810.win 

06/25/98 

Gap9811.win 

06/25/98 

Gap9812.win 

06/26/98 

Gap9813.win 

06/26/98 

Gap9814.win 

06/26/98 

Gap9815.win 

06/27/98 

Gap9816.win 

06/27/98 

Gap9817.win 

06/27/98 

Launch  time  LT 

Max  Altitude 

2215 

31.08 

0025 

32.35 

2332 

30.07 

2130 

30.09 

2320 

28.12 

1548 

30.27 

2140 

19.50 

2345 

23.44 

2130 

27.63 

2359 

30.03 

1935 

30.02 

2133 

9.43 

2241 

30.18 

1730 

22.24 

1926 

29.79 

2130 

30.27 
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APPENDIX  2 

Data  Additional  information 

Question  from  Y.  Louvet  on  19th  June  2001  and  answer  from  Dr.  Jumper 

>As  you  know  it  is  always  very  difficult  to  interpret  measurements  when 
>haven’t  to  take  part  in  the  trial  measurements,  so  we  have  again  a  lot  of 
>questions  from  Dr.  Dalaudier. 

> 

>1-  The  data  are  difficult  to  correlate  because  we  do  not  have  the  exact 
>dating  of  the  measurement.  Have  you  kept  the  dating  of  measurements  and 
>could  you  send  it  to  us? 

>We  have  all  dates  and  times  of  launch.  The  balloon  goes  up  at  5  to  7  m/s, 

>so  there  is  a  lag  as  it  progresses.  I  have  attached  the  time  sheet  summary 
>below  (GAP98T.pdf).  The  time  and  date  are  "Local"  Gap  time,  NOT  Universal 
>Time.  If  you  have  a  problem  with  the  Adobe  PDF  document,  let  me  know  and  I 
>will  send  you  a  text  copy.  (In  fact,  I  noticed  the  math  fonts  are  not  very 
>good  -  if  anyone  wants  a  good  copy,  I  can  send  it  by  mail.) 

>2-  In  addition  we  have  some  difficulties  in  correlation  of  the  altitude 
>data.  How  is  calculated  the  data  of  altitude  from  the  file  *.win?  is  it 
>altitude  coming  from  data  GPS  or  is  it  from  calculated  altitude  by 
integration  of  vertical  winds  measurements  of  the  Vaisala  probe? 

>Altitudes  are  computed  the  standard  way  for  a  radiosonde.  The 
"hypsometric" 

>equation  is  used,  that  is,  the  hydrostatic  equation  with  humidity 
considered 

>in  the  equation.  The  radiosonde  measures  pressure,  temperature,  and 
>humidity  (and  we  get  wind  speed). 

>3-  The  probe  Vaisala  gives  a  frequency  of  measurements  which  is  not  usual 
>acts  of  a  probe  RS90? 

> 

>  4-  In  the  *.text  file,  how  is  measured  or  calculated  the  pressure  ? 

>Pressure  is  measured.  Altitude  is  computed. 

>5-  It  is  always  interesting  to  know  the  processing  of  the  signal  carried 
>out  of  differential  measurement  with  the  electronics  bandwidths.  How  the 
differential  temperature  signal  is  it  treated  and  how  much  is  the 
bandwidth 

>of  electronic  used?. 

>The  attached  paper  "  AIAA  Sources  of  Error"  below  discusses  the  instrument 
>and  the  measurements  along  with  references.  (It  also  discusses  possible 
>sources  of  error).  The  differential  temperature  is  measured  with  some  very 
>fine  wire  probes,  3.45  micro  meters  in  diameter.  The  time  constant  of  the 
>wires  are  nominally  above  250  hz  at  maximum  altitude  -  faster  on  the 
ground. 

>The  band  pass  of  the  filters  is  between  1000Hz  and  0.5Hz.  The  amplified 
and 

>filtered  signal  is  then  passed  through  an  RMS  chip  with  a  nominal  time 
>constant  of  3.75s.  As  discussed  in  the  paper,  this  was  found  to  actually 
depend  on  whether  the  signal  was  rising  or  falling.  When  falling,  the  time 
>constant  increases  to  7.5s. 

> 
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Question  from  Dr.  Louvet  on  12th  July  2001 
Answer  from  Dr.  Jumper  on  27th  July 

1)  Firstly  about  very  interesting  publication  AIAA  Sources  of 
Error.pdf  unfortunately  as  you  be  afraid  of  a  numerous  problems  happened 
on  math  fonts.  Please  can  you  send  a  text  copy,  see  below  address.  Thank 

. See  attached  .pdf  file.  Hopefully  the  equations  will  be  readable  this 

time... 

As  you  know  atmospheric  stability  study  requires  very  good  data 
correlation  between  measurements.  So  we  try  to  understand  time  acquisition, 
pressure  measurement  and  calculated  hydrostatic  altitude  correlation.  We 
are  a  little  bit  surprise  by  data  can  you  help  us. 

2)  The  ratio  of  pressure  variation  by  pressure  seems  excellent  and  constant 
as  altitude.  What  is  the  sensor  reference? 

What  is  VAISAL  probe  reference  and  work  mode  acquisition  choosing 

...  It  is  the  Vaisala  pressure  sensorfor  the  RS80.  We  will  FAX  the 
specifications  to  you  along  with  the  specifications  on  the  Ozone  sensor  (FAX 
to  Laurent)... 

3)  The  dat  files  show  that  time  acquisition  is  not  constant  and  is  very 
different  as  time  difference  products  from  cutting  time,  (for  exemple  if 
time  is  only  measured  in  seconds  and  if  acquisition  frequency  is  1 ,2s  we 
have  time  data  recorded  steps  of  1  and  2  seconds). 

...  For  some  strange  reason,  Vaisala  rounds  the  times  in  the  files  to  the 
nearest  second,  even  though  it  shows  several  decimal  places.  Those  times 
can  be  used  as  is,  accepting  a  +/-  .5  second  error,  or  the  times  can  be 
replaced  by  multiples  of  the  sample  interval.  If  you  do  the  latter,  you 
have  to  be  careful,  because  sometimes  a  line  is  lost  in  the  transmission 


In  txt  file  it  seems 

that  measurements  are  taken  one  or  two  or  three  time  of  fundamental 
acquisition  frequency  as  we  see  in  *.dat  file.  Please  do  you  know  why  ? 

...  In  the  .txt  file,  we  present  data  for  approximately  equal  values  of 
time,  and  interpolate  if  required.  There  is  no  time  given,  only  altitude  in 
km ... 

4)  In  *.TXT  we  have  several  values  of  Ozone  pressures  (P303,  P403, ..)  can 
you  remenber  me  what  is  the  difference  and  where  come  from  these  different 
values? 

...The  radiosonde  actually  has  5  spare  channels  that  we  can  use  for  the 
additional  data,  and  a  couple  of  reference  voltages  .  All  the  03  data  comes 
from  the  same  instrument,  it  is  just  sampled  5  times  if  it  is  the  only 
instrument.  If  we  have  the  thermosonde  also  attached  to  2  or  3  of  the  spare 
channels,  we  show  2  or  3  03  data  points  and  3  or  2  thermosonde  data  points. 
Typically  we  then  use  only  one  of  the  channels  for  analysis  ... 
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5)  About  wind  angle  we  have  supposed  it  is  right  to  meteo  convention  which 
gives  angles  where  the  wind  come  from  and  not  the  angle  of  wind  vector 
goes. 

...You  are  correct,  the  meteo  convention  is  used,  it  is  the  angle  the  wind 
is  coming  from.  0  is  from  the  North,  90  is  from  the  East.... 

6)  About  time.  In  your  mail  of  the  19  th  June  you  have  confirmed  according 
to  GAP98T.pdf  file  that  the  time  is  in  Local  Time.  Unfortunetly  in 
dataform.doc  file  send  with  txt  and  win  file  it  is  write  GMT  time  and  it 

is  the  same  values  ?  Can  you  confirm  the  right  time? 

...  In  the  table  of  all  launches,  we  give  the  local  date  and  time  at 
launch.  It  is  local  time  that  is  sometimes  shown  on  graphs.  There  is  no 
time  in  the  .txt  and  .win  files  -  the  first  column  is  altitude  (km).  The 
.met  file  does  have  time  to  the  nearest  second,  with  an  arbitrary  zero  time 
near  the  actual  launch  time.  We  sent  the  .dat  files  so  that  you  could  tell 
what  time  (seconds  after  launch)  was  associated  with  each  altitude.  You  can 
add  those  seconds  to  launch  time  to  get  actual  time.  Remember  any  launch 
time  we  provide  is  local.  I  do  not  know  why  someone  said  there  was  GMT  in 
the  .txt  and  .win  files.  There  is  no  time  at  all  in  those  files  ... 

7)  If  possible  we  are  very  interested  of  all  data  from  the  descent .  We 
know  that  at  the  beginning  descent  the  data  are  certainly  no  right,  but 
have  you  any  think  ? 

...  We  do  not  know  much  about  the  descent  of  that  one  flight.  The  ozone 
concentrations  agree  fairly  well,  so  we  showed  them  in  the  paper... 

8)  Have  you  the  geographic  path  as  a  function  of  time  for  each  launch.  Can 
you  send  it. 

...  We  did  not  compute  it.  We  sometimes  assume  that  we  are  totally  in 
equilibrium  with  the  horizontal  wind  and  compute  the  position  from  the  wind 
file,  but  we  did  not  accomplish  that  for  these  flights  ... 

9) Good  new.  Preliminary  analysis  of  ozone  gradient  shows  a  very  good 
correlation  between  potential  temperature.  It  means  that  spatial  resolution 
of  03  sensor  was  better  than  100  m  in  the  altitude  range  of  18  to  22  km 
high,  have  you  some  idea  about  time  response  of  o3  probe  as  function  of 
altitude. 

...  I  do  not  remember  the  time  response  as  a  function  of  altitude.  I 
thought  that  it  was  associated  with  the  volume  of  the  pump,  which  does  not 
change,  but  I  might  be  wrong.  We  will  look  into  this  and  get  back  to  you  if 
we  find  anything.... 
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APPENDIX  3 
FIRST  SET  OF  GRAPHS 
RUNS  9,10, 11, 12, 14, 15  AND  16 


For  each  run  we  plotted  9  figures,  in  the  following  order: 


Figure  1  :  Wind  components  UU  (West  -East )  et  VV  (South-North)  in  the  [0-35]  km  altitude  range. 

Figure  2  :  Zoom  in  the  [1 5-1 8]  km,  [18-22]  km  and  [22-30]  km  altitude  ranges. 

Figure  3  :  Plots  of  UU,  VV,  measured  (T)  and  potential  (Tpotent)  temperatures,  ozone  concentration  in 

the  [0-35]  km  altitude  range. 

Figure  3bis  :  Same  as  3  in  the  [15-18]  km,  [18-22]  km  and  [22-30]  km  altitude  ranges. 

Figure  4  :  Wind  components  UU  (W-E)  et  VV  (S-N);  03  concentration  and  gradient  and  N2  (square 

Brunt-Vaisala  pulsation  multiplied  by  1 0000)  in  the  [1 8-22]  km  altitude  range. 

Figure  5  :  Same  as  Figure  4  but  on  two  graphs. 

Figure  6  :  Curves  of  P  (pressure),  T,  Tpotent,  03  and  %humidity  in  the  [0-35]  km  altitude  range. 

Figure  7  :  Graphs  of  shear  (wind  gradient),  N2  et  cn2  (runs  9,  1 0,  1 1 ,  and  14)  in  the  [0-35]  km  altitude 

range. 

Figure  8  :  Same  as  Figure  7  but  shear,  N2  et  cn2  (runs  9, 1 0, 1 1 ,  and  14)  in  the  [18-22]  km  altitude  range 
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Gap98-RUN  09  Fig3  24/08/01^ 
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Wind  Gap98  RUN  11  Fig  2  24/08/01 YL 
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SECOND  SET  OF  GRAPHS 
RUNS  9, 10, 11, 14 


We  have  nine  figures  for  each  run. 

Figures  X.4.1,  X.4.2  and  X.4.3  represent  ozone  gradient  and  concentration,  wind  components  U  and 
V  and  N2  (noted  L2).  Attention:  L2  is  in  (rad/s)2  *1000. 

Figures  X.8.1,  X.8.2  and  X.8.3  show  the  shear  (wind  gradient  in  [m/s]/km  and  divided  by  10),  cn2 
and  Richardson  number. 

Figure  9  shows  cn2  as  a  function  of  shear. 


We  have  detected  some  missing  link  on  wind  data 
(RUN  9  from  15  to  17  km 
RUN  10  none, 

RUN  1 1  from  3  to  7  km,  from  18.5  to  19.3  and  19.7  km  and  27  km 
RUN  14  at  19.3  km  and  18.7  km 


For  each  run  we  plotted  9  figures  in  the  following  order: 

Figure  1  :  Wind  components  UU  (West  -East )  et  VY  (South-North)  in  the  [0-35]  km  altitude  range. 

Figure  2  :  Zoom  in  the  [15-18]  km,  [18-22]  km  and  [22-30]  km  altitude  ranges. 

Figure  3  :  Plots  of  UU,  VV,  measured  (T)  and  potential  (Tpotent)  temperatures,  ozone  concentration  in 

the  [0-35]  km  altitude  range. 

Figure  3bis  :  Same  as  3  in  the  [15-18]  km,  [18-22]  km  and  [22-30]  km  altitude  ranges. 

Figure  4  :  Wind  components  UU  (W-E)  et  VV  (S-N);  03  concentration  and  gradient  and  N2  (square 

Brunt-Vaisala  pulsation  multiplied  by  10000)  in  the  [18-22]  km  altitude  range. 

Figure  5  :  Same  as  Figure  4  but  on  two  graphs. 

Figure  6  :  Curves  of  P  (pressure),  T,  Tpotent,  03  and  %humidity  in  the  [0-35]  km  altitude  range. 

Figure  7  :  Graphs  of  shear  (wind  gradient),  N2  et  cn2  (runs  9,  10,  11,  and  14)  in  the  [0-35]  km  altitude 

range. 

Figure  8  :  Same  as  Figure  7  but  shear,  N2  et  cn2  (runs  9,  10,  1 1 ,  and  14)  in  the  [18-22]  km  altitude  range. 
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DESCRIPTION  OF  MEASUREMENTS  FROM  THE  SERVICE  D’AERONOMIE 

SAOZ 

The  SAOZ  instrument  [4-5]  consists  of  a  UV-visible  spectrometer  allowing  measurement  of  the 
vertical  profiles  of  03,  N02,  OCIO,  BrO  and  10  by  solar  screening  during  the  rise  of  the  balloon  in  twilight 
from  a  ceiling  close  to  30  km.  Controlled  by  a  computer,  it  includes  a  pressure  and  temperature 
measurement,  localization  by  Global  Positioning  System  (GPS)  and  an  Argos  transmitter  for  its  recovery  on 
the  ground  after  parachute  descent.  The  total  weight  of  20  kg  is  designed  for  frequent  and  cheap  flights 
under  5,000  and  10,000  m3  balloons  developed  by  the  National  Center  of  Space  Studies  (CNES).  Numerous 
flights  have  been  carried  out  since  1991  in  France,  Spain,  Norway,  Sweden  and  Brazil  to  study  the  ozone 
destruction  and  the  chemistry  of  chlorine,  bromine  and  iodine. 

The  list  of  SAOZ  flights  is  given  on  web  site  http://www.aerov.jussieu.fr/~fgoutail/FlightList.html. 
Two  flights  including  ozone  measurements  were  made  in  June  1998  from  the  Gap  launch  base,  on  June  20  at 
sunrise  and  June  23  at  sunset  (see  Figure  2).  Dr.  Goutail  of  the  Service  d’Aeronomie  is  in  charge  of  the 
instrument. 

EEC  Ozone  Probe 

The  Observatoire  de  Haute  Provence  station  makes  an  atmospheric  ozone  survey  by  EEC  probe 
every  Wednesday.  A  survey  was  made  on  Wednesday  June  24  during  the  trial  runs.  These  measurements 
were  provided  by  Mr.  Vialle  who  is  the  measurement  PI  at  the  Service  d’Aeronomie. 

Stratospheric  O3  Lidar 

The  stratospheric  ozone  lidar  used  at  OHP  is  a  UV  Dial  Lidar  [6-7]  emitting  at  two  wavelengths,  one 
of  which  is  absorbed  by  ozone.  The  laser  pulses  are  scattered  by  atmospheric  particles  at  various  altitudes. 
Measurements  are  made  by  time-of-flight  detection  of  the  radiation  scattered  by  the  reflecting  layers. 
Comparison  of  the  range  signals  from  the  two-wavelength  scattering  laser  gave  an  ozone  measurement 
integrated  on  the  probed  column. 

Measurements  were  made  on  the  nights  of  June  19  and  between  June  21  and  June  28. 

The  validity  range  of  measurements  in  the  10  km  to  45  km  altitude  range  is  supplied  in  the  data  file 
header  and  must  be  taken  into  account.  The  person  in  charge  of  these  measurements  is  Dr.  Godin  of  the 
Service  d’Aeronomie. 

Rayleigh  Lidar 

A  Rayleigh  lidar  [8-12]  provided  temperature  measurements  with  a  spatial  resolution  of  75  m  in  the 
altitude  range  from  25  km  to  70  km.  The  system  detects  scattering  of  laser  pulses  by  Rayleigh  molecules. 
The  time-of-flight  of  the  backscattered  signal  gives  the  distance  from  the  probed  layers.  The  measurements 
give  access  to  atmospheric  density  which  can  be  converted  to  temperature  using  the  hydrostatic  coefficients. 

Measurements  below  30  km  can  be  wrong  because  of  aerosols.  Measurements  were  made  nightly  on 
June  19  and  from  June  21  to  June  29.  The  measurement  times  are  given  in  Figure  3.  The  measurements 
lasted  approximately  5  hours,  beginning  at  sunrise  Dr.  Hauchecorne  of  the  Service  d’Aeronomie  was  the 
leader  of  these  measurements. 
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50  years  at  the  cutting  edge  of  civil  and  military 
aerospace  research 

ONERA  has  been  involved  in  the  development  of  a  number  of  landmark 
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reputation  for  scientific  and  technical  excellence. 

A  revamped  onera  steps  up  the  pace 

ONERA  has  revamped  its  structures  to  keep  pace  with  the  fast-changing 
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